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Abstract 

Variational Monte Carlo wave functions, obtained from a realistic Hamilto- 
nian consisting of the Argonne v\% two-nucleon and Urbana-IX three-nucleon 
interactions, are used to calculate the 6 Li ground-state longitudinal and trans- 
verse form factors as well as transition form factors to the first four excited 
states. The charge and current operators include one- and two-body com- 
ponents, leading terms of which are constructed consistently with the two- 
nucleon interaction. The calculated form factors and radiative widths are in 
good agreement with available experimental data. 

PACS numbers: 21.45. +v,25.30.Bf,25.30.Dh,27.20.+n 
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Calculations of the 6 Li elastic and inelastic form factors have relied in the past on rel- 
atively simple shell-model or a-d 0] cluster wave functions. These calculations have 
typically failed to provide a satisfactory, quantitative description of all measured form fac- 
tors. More phenomenologically successful models have been based on aNN J5|-[7| clusteriza- 
tion, on extensions of the basic spherical-cluster a-d model in which the deuteron is allowed 
to deform |J, or on large-space multi-Ziu; shell- model approaches 0. However, while these 
models do provide useful insights into the structure of the A=6 nuclei, their connection with 
the underlying two- (and three-) nucleon dynamics is rather tenuous. 

The 6 Li form factor calculations we report on here are within the context of a realistic 
approach to nuclear dynamics based on two- and three-nucleon interactions — the Argonne 
t>i8 [[0J and Urbana IX fTl] interactions, respectively, or AV18/UIX model — and consistent 
two-body charge and current operators p2|~P^]. Up until very recently, calculations of this 
type were limited to the A=2-4 systems, as reviewed in Ref. |15|1 . Indeed, the deuteron 
structure functions and threshold electrodisintegration, the trinucleon charge and magnetic 
form factors, and a charge form factor have been the observables of choice for testing the 
quality of interactions and associated two-body currents. However, the availability of real- 
istic six-body wave functions for the ground and low-lying excited states of 6 Li |1| makes 
it now possible to extend and test our understanding of the electromagnetic structure of 
nuclei in a new regime — that of p-shell nuclei — and to verify to what extent the inability 
of reproducing simultaneously the observed elastic and transition form factors is due to the 
inadequacy of cluster or shell-model wave functions. 

The AV18/UIX model reproduces the experimental binding energies and charge radii 
of 3 H, 3 He, and 4 He in numerically exact calculations, based on the Faddeev W7\, pair- 

and Green's function Monte Carlo (GFMC) HI 
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correlated hyperspherical harmonics 
methods. For A=6 systems the GFMC results are somewhat underbound compared to the 
experimental ground states, by 2% in 6 Li and 5% in 6 He and 6 Be. The best variational 
Monte Carlo (VMC) energies are an additional 10% above the GFMC results. However, the 
known excitation spectra are well reproduced by both VMC and GFMC calculations. These 
include, in order of excitation, the states with spin, parity, and isospin assignments, ( J 71 "; T), 
of(3+;0), (0+;l), (2+;0), and (2+; 1) 0. 

The variational wave function for A=6 nuclei used here is the trial wave function, ^t, 
that serves as the starting point for the GFMC calculations. It has the general form 



i + E 
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where Uy and U™ 1 are two- and three-body correlation operators and the Jastrow wave 
function is given by 
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The S and A are symmetrization and antisymmetrization operators, respectively. The 
central pair and triplet correlations fxy(j"ij) and f^ k are functions of relative position only; 
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the subscripts xy denote whether the particles are in the s- or p-shell. The |$ 6 (L5 J MTT 3 )) 
is a single-particle wave function with orbital angular momentum L and spin S coupled to 
total angular momentum J, projection M, isospin T, and charge state T 3 : 

\$ 6 (lsjmtt 3 ) 12M:56 ) = l$4(oooo) 1234 </>f (R a5 )<i>i s (Ra e ) 

{[Y lmi (n a5 )Y lm ,(n a6 )] LML x l X5 (±m s )xe(lrn' s )] SMs } jm 
xH{\^X¥z)\tt,) ■ (3) 

Particles 1-4 are placed in an a core with only spin-isospin degrees of freedom, denoted by 
$4(0000), while particles 5-6 are placed in p-wave orbitals <f)p S (R a k) that are functions of 
the distance between the center of mass of the a core and particle k. Different amplitudes 
Pls are mixed to obtain an optimal wave function; for the (l + ;0) ground state of 6 Li we 
mix /3 01 , Piq, and (3 2 i terms, while both the (3 + ;0) and (2 + ; 0) states are "stretch "states 
and use only /3 2 %. For the (0 + ; 1) excited state the /3 00 and 0u amplitudes contribute, and 
for (2 + ; 1) the wave function is constructed from (3 20 and 0n terms. 
The two-body correlation operator EL,- is defined as: 



u* = E 

p=2,6 



II fijk( T ik, r jk) 

k^i,j 
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where the 0\~ ' = • Tj, cr^ ■ crj, a \ ■ (TjT , • Tj, S^, and SyTj • t 3 -. The six radial functions 
f ss {r) and u p=2 fi (r) are obtained from approximate two-body Euler-Lagrange equations with 
variational parameters []I9]. The / sp and /^f correlations are similar to / ss for small sepa- 
rations, but include long-range tails. The parameters used in constructing these two-body 
correlations, as well as the description of the three-body correlation operator U™ 1 and the 
operator-independent three-body correlations fh k and ffj k are given in Ref. [Tj| . 



Energy expectation values are evaluated using a Metropolis Monte Carlo algorithm [|19 
The VMC results for the ground and low-lying excited states of 6 Li are compared to the 



GFMC and experimental energies |20[ in Table I. The ground state is underbound by nearly 
5 MeV compared to experiment, and is only 0.1 MeV more bound than the corresponding 
4 He calculation (26.9 MeV). This is above the threshold for breakup of 6 Li into an a and 
a deuteron; in principle, it should be possible to lower the variational energy at least to 
that threshold, but the wave function would be greatly spread out. We have chosen to 
constrain our parameter search to keep the rms point-nucleon radius for the ground state 
near the experimental value of 2.43 fm. Despite the large energy deficit compared to the 
GFMC calculation, the VMC and GFMC one-body densities in 6 Li are virtually identical. 
However, the two-body GFMC densities are somewhat larger near their peak at pa 1 



fm [16 



The nuclear charge and current operators consist of one- and two-body terms. We here 
summarize their most important features, and refer the reader to Refs. []15| , pT] for a listing 
of the explicit expressions. The two-body current operator has "model- independent "and 
"model-dependent "components, in the classification scheme of Riska [Q. The model- 
independent terms are obtained from the charge-independent part of the AV18, and by 



construction [123 j satisfy current conservation with this interaction. The leading operator 



is the isovector "7r-like "current obtained from the isospin-dependent spin-spin and tensor 
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interactions. The latter also generate an isovector "p-like "current, while additional model- 
independent isoscalar and isovector currents arise from the isospin-independent and isospin- 
dependent central and momentum-dependent interactions. These currents are short-ranged 
and numerically far less important than the 7r-like current. 

The model-dependent currents are purely transverse and therefore cannot be directly 
linked to the underlying two-nucleon interaction. The present calculation includes the 
isoscalar pivy and isovector uTf-f transition currents as well as the isovector current associated 
with excitation of intermediate A-isobar resonances. The prc-f and uTf-f couplings are known 
from the measured widths of the radiative decays p — > 777 [|23J and uj — ► 777 respectively, 
while their momentum-transfer dependence is modeled using vector-meson-dominance. The 



Ml jNA coupling is obtained from an analysis of data in the A-resonance region p6 



Monopole form factors are introduced at the meson-baryon vertices with cutoff values of 
A 7r =3.8 fin" 1 at the vriViV and vriVA vertices and A p =A w =6.3 mi" 1 at the pNN and uNN 
vertices. 

While the main parts of the two-body currents are linked to the form of the two-nucleon 
interaction through the continuity equation, the most important two-body charge operators 
are model-dependent, and should be considered as relativistic corrections. Indeed, a consis- 
tent calculation of two-body charge effects in nuclei would require the inclusion of relativistic 
effects in both the interaction models and nuclear wave functions. Such a program is just at 
its inception for systems with A > 3. There are nevertheless rather clear indications for the 
relevance of two-body charge operators from the failure of the impulse approximation (IA) 



in predicting the charge form factors of the three- and four-nucleon systems [L5||. The model 



commonly used [|13|,|nj includes the it-, p-, and w-meson exchange charge operators with 
both isoscalar and isovector components, as well as the (isoscalar) piry and (isovector) umj 
charge transition couplings, in addition to the single-nucleon Darwin-Foldy and spin-orbit 
relativistic corrections. It should be emphasized, however, that for q<5 fm -1 the contribu- 
tion due to the 7i-exchange charge operator is typically an order of magnitude larger than 
that of any of the remaining two-body mechanisms and one-body relativistic corrections. 

We have calculated longitudinal F^q) and transverse Fj<(q) form factors of the 6 Li ground 
state as well as transitions from this to the first four excited states. The Coulomb (CJ) 
multipoles contributing to Fi(q) are obtained from matrix elements of the charge operator, 
p(q), while the electric (EJ) and/or magnetic (MJ) multipoles contributing to Fx(q) are 
obtained from matrix elements of the current operator, j(q), using standard formulas p7| . 

Our elastic form factors F^q) and Ff{q) are compared with the experimental val- 
ues |8|,p^-|30| in Fig- 0- Since the 6 Li ground state is (l + ;0), both CO and C2 multipoles 
contribute to F L (q), while only the Ml operator contributes to F T (q). The results obtained 
in both IA and with inclusion of two-body corrections in the charge and current operators 
(IA+MEC) are displayed, along with the statistical errors associated with the Monte Carlo 
integrations. The Fi(q) is in excellent agreement with experiment; in particular, the two- 
body contributions (predominantly due to the 7r-like charge operator) shift the minimum to 
lower values of q, consistent with what has been found for the charge form factors of the 
hydrogen and helium isotopes | 15| . The C2 contribution is much smaller than CO below 
3 fm -1 , as shown in Fig. [I], and at low q is proportional to the ground state quadrupole 
moment. Our prediction for the latter is —0.23(9) fm 2 , larger (though with a 50% statistical 
error) in absolute value than the measured value of —0.08 fm 2 , but with the correct (nega- 
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tive) sign. Cluster models of the 6 Li ground state generally give large, positive values for the 
quadrupole moment, presumably due to the lack of D- waves in the a, and the consequent 
absence of destructive interference between these and the D-wave in the a-d relative mo- 
tion [J7|]. For q > 3 fm -1 , the C2 contribution becomes dominant, and the shoulder seen in 
the data is entirely due to this component, which has frequently been omitted from cluster 
models. 

The experimental -Pr(g) is well reproduced by our calculations in the first peak at q — 
0.5 fm -1 , but the zero comes a little too early and the second peak at q = 2 fm -1 is 
somewhat overpredicted. Since the 6 Li ground state has T=0, only isoscalar two-body 
currents contribute to Fxlq); the associated contributions are small at low q, but increase 
with q, becoming significant for q > 3 fm -1 , beyond the range of present data. The calculated 
magnetic moment is 0.829 in IA and 0.832 hn with two-body currents, about 1% larger 
than the experimental value, which is close to that of a free deuteron. 

The calculated longitudinal inelastic form factor to the (3 + ;0) state is found to be in 
excellent agreement with experiment |2^j31|j3^] , as shown in Fig. |2|. This transition is induced 
by C2 and C4 operators, and thus the associated form factor Fi(q) behaves as q 4 at low 
q. The two-body contributions only become important for q > 2 fm -1 , but do improve the 
agreement with data. We also show our prediction for the much smaller transverse form 
factor F T (q). The calculated radiative width of the (3 + ;0) state is 3.38(9) x!0~ 4 eV in both 



IA and with MEC, compared to the experimental value of (4.40 ± 0.34) x 10 4 eV [|33 



Good agreement is also found with the experimental values [0,^2j for the transverse 
inelastic form factor to the state (0 + ;l), as shown in Fig. 0. This is an isovector magnetic 
dipole transition and, as expected, is significantly influenced, even at low values of q, by 
two-body contributions, predominantly by those due to the 7r-like current operator. The 
predicted radiative width is 7.49(2) eV in IA and 9.06(7) eV including MEC, compared with 
the experimental value 8.19 ± 0.17 eV. Thus the isovector two-body current contributions 
increase the 7-width by 20%. 

The calculated longitudinal and transverse inelastic form factors to the (2 + ;0) state are 
also shown in Fig. [2]. The contributing multipole operators are C2 for the longitudinal 
transition, and Ml, E2, and M3 for the transverse transition. The F^q) is comparable in 
magnitude to that for the (3 + ;0) state, but has not been measured to date; the Fxiq) is 
again much smaller. The corresponding 7-width is calculated at 8.0(5) xl0~ 3 eV both in IA 
and with MEC. This result just overlaps the experimental value of (5.4 ± 2.8) xlO -3 eV |33 



Finally, we show the inelastic form factors to the (2 + ;l) state. For this isovector tran- 
sition, the F T (q) is much larger than the F L (q). The experimental data [fj| |3~2^|33|l are well 
reproduced by the calculation, with the two-body currents contributing significantly at all 
q values. However the calculated 7-width of 0.050(9) eV in IA and 0.075(26) eV with MEC 
is several times smaller than the reported experimental value of 0.27 ± 0.05 eV. 

To summarize, we have presented the first ah initio microscopic calculations of 6 Li elastic 
and transition form factors, based on six-body VMC wave functions obtained from realistic 
interactions and a consistent, realistic nuclear electromagnetic current operator. We do 
not expect that use of the more accurate GFMC wave functions will lead to significantly 
different predictions, since the VMC and GFMC one- and two-body densities have been 
found to be quite close ]16j. Inclusion of the contributions from two-body charge and current 



operators brings theory into significantly better agreement with the experimental data. Thus 
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nuclear many-body theory appears to provide a quantitatively satisfactory description of 
the electromagnetic structure of both s- and p-shell nuclei for a wide range of momentum 
transfers. 
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TABLES 







TABLE I Bindine 
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FIGURES 



FIG. 1. Calculated longitudinal and transverse elastic form factors of the 6 Li ground state 
are shown in impulse approximation (IA) and with two-body charge and current operators added 
(IA+MEC) as filled symbols with Monte Carlo statistical error bars. The Coulomb monopole (CO) 
and quadrupole (C2) contributions to the longitudinal form factor are also shown by the dashed 
(IA) and solid (IA+MEC) lines. Data are from Refs. gpHHO]. 



FIG. 2. Calculated longitudinal and transverse transition form factors to the first four (J 71 ";! 1 ) 
excited states of 6 Li are shown in impulse approximation (IA) and with two-body charge and 
current operators added (IA+MEC). The largest form factor in each case is shown as a point 
with its statistical error bars; the smaller form factor (if any) is shown by dashed (I A) and solid 
(IA+MEC) lines. Data are from Refs. [||]J^j3l]-|3|. 



9 



Fig.l Wiringa & Schiavilla 




Fig. 2 Wiringa & Schiavilla 
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